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Abstract. A (p, <j)-ary chain is a special type of chain partition of integers 
with parts of the form p a q b for some fixed integers p and q. In this note, we 
are interested in the maximal weight of such partitions when their parts are 
distinct and cannot exceed a given bound m. Characterizing the cases where 
the greedy choice fails, we prove that this maximal weight is, as a function of 
m, asymptotically independent of max(p, q), and we show how to compute its 
value in logarithmic time. 



1. Introduction 

Let p, q be two fixed integers, and let E — {p a q b : (a, b) e N 2 } be endowed with 
the divisibility order, i.e. x >z y y \ x. A (p, q)-ary chain is a finite non- 

increasing sequence in E. For example, (72, 12, 4, 4, 1) is a (2, 3)-ary chain, whereas 
(72, 12, 4, 3, 1) is not since 4^3. We define the weight of a (p, <?)-ary chain as the 
sum of its terms, i.e. an expression of the form 

w = ^p Q, g bl , where p a 'q b ' h p a '+^q b w for i > 1. (1) 
i>i 

Expansions of this type have been proposed and successfully used by Dimitrov in 
the context of digital signal processing or cryptography under the name double-base 
number system. For more details, see [TJ [2] and the references therein. 

From a different point of view, a (p, <?)-ary chain can be seen as a partition of its 
weight, where the parts are restricted to the set E and constrained by a divisibility 
condition. Surprisingly, works on integer partitions with divisibility constraints 
on the parts are very scarce. Erdos and Loxton [3] considered two types of such 
unconventional partitions, called chain and umbrella partitions, and obtained "some 
rather weak estimates for various partition functions". More recently, motivated by 
some theoretical questions behind Dimitrov's number system, the second and third 
authors refined some of Erdos and Loxton's earlier results in a paper [4] dedicated 
to strictly chained (p,q)-ary partitions. A strictly chained (p, q)-axy partition, or 
(p, q)-SCP for short, is a decreasing (p, (?)-ary chain, i.e. it has distinct parts. In 
this note, we are interested in the maximal weight of a (p, g)-SCP whose parts are 
bounded by some given integer m. In particular, assuming p < q, we prove that 
this maximal weight asymptotically grows as pmj (p — 1), independently of q. 

The heaviest (p, g)-SCP whose first part is given may be computed using a greedy 
strategy: successively take the next greatest part satisfying the divisibility condi- 
tion. Nevertheless, given a bound m > on the parts, determining how to best 
select the first part is not immediate and the greedy approach appears to fail as a 
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general rule. That is, choosing the largest part less than or equal to m does not 
always provide a partition of maximal weight. These facts are established in Sec- 
tions and [3] among other preliminary definitions, examples, and results. The cases 
where the greedy choice fails are fully characterized in Section 01 while Section O is 
devoted to the asymptotic behavior of the maximal weight as a function of bound 
to. Finally, it is explained in Section how to compute a best choice for the first 
part, thus the maximal weight, in O(loglogm) time. 

It is worthwhile to point out that our results still hold for any pair of multiplica- 
tively independent algebraic numbers p, q. 

2. Preliminaries 

Let to be a positive integer, and let G(to) denote the maximal weight of a 
(p, g)-SCP whose greatest part does not exceed to. For example, with p = 2 and 
q = 3, the first values of G are: 1, 3, 4, 7, 7, 10, 10, 15, 15, 15, 15, 22, 22, 22, 22, 
31, 31, ... 

In the following, we shall assume w.l.o.g. that p < q. Notice that the case p = 1 
is irrelevant since G(to) is simply the sum of all the powers of q less than or equal 
to to. More generally, and for the same reason, we shall consider that p and q are 
not powers of the same integer, that is, they are multiplicatively independent (see, 
e.g., [5] Th. 2.5.7). Under this assumption, the first values of G(m) may be quickly 
computed with the help of the following formula. 

Proposition 1. For to 6 N*, let G(m) denote the largest integer that can be 
expressed as a strictly chained (p,q)-ary partition with all parts less than or equal 
to to. Assume that G(to) =0 if to N. Then, we have G(l) = 1, and for to > 1 

G(to) =max(G(m.- + pG(m/p),l + qG(m/q)). (2) 

Proof. Let A be a partition of weight G(to) whose parts are all less than or equal 
to to. First, notice that A must contain part 1 by definition of G(m). If to is not an 
element of E, then G(to) — G(to — 1). Otherwise, it suffices to note that removing 
part 1 in A produces a partition whose parts are all divisible by either p or q. □ 

Computing G(to) with relation ^ requires O(logTO) steps in the worst case: 
simply note that, for all to, in at most p — 1 baby-steps, i.e. G(to) = G(to — 1), 
one gets an integer that is divisible by p. Formula ^ may also be adapted to 
provide a (p, q)-SCP of weight G(to). Nevertheless, it does not give any idea about 
the asymptotic behavior of G. Moreover, we shall see in Section [5] how to compute 
G(to) and a (p, g)-SCP of weight G(m) in O(loglogm) time. 

A natural graphic representation for (p, g)-SCPs obtains by mapping each part 
p a q b G E to the pair (a, b) G N 2 . Indeed, with the above assumptions on p and q, 
the mapping (a, b) i— > p a q b is one-to-one. Since the parts of a (p, g)-SCP are pairwise 
distinct by definition, this graphic representation takes the form of an increasing 
path in N 2 endowed with the usual product order. This is illustrated in Figure [1] 
with the ten (2, 3)-SCPs containing exactly six parts and whose greatest part equals 
72 = 2 3 3 2 . Note that a (p, q)-SCP with largest part p a q b has at most a + b+1 parts, 
and there are exactly ( a ^ h ) of them with a maximum number of parts. 

With this representation in mind, one is easily convinced that the heaviest 
(p, q)-SCP with first part p a q b looks like the top left (p, g)-SCP in Figured) This is 
formalized in the following lemma. 
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Figure 1. The set of (2,3)-SCPs with 6 parts and whose largest 
part equals 2 3 3 2 = 72. 



Lemma 1. Given a, b G N, the heaviest (p,q)-SCP with first part p a q b is the one 
whose parts are the elements of the set {q 1 : < i < b) U {q b p l '■ < i < a}. 



Proof. Consider a (p, q)-SCP A = (\)i =1 with greatest part A x = p a q b . Let A^ 



p^q 



b, 



If a, + hi > b then define \' t = p ^+^-b q b^ otherwise let A^ = q a > 



-bi 



Note 



that X[ = p a q again. Since sequence (ai + bi)^ =l is decreasing, (A^)f =1 is also a 
(p, q)-SCP. Since p < q we have A^ > A^ for all i, with equality if, and only if, the 
parts in A form a subset of {q l : < i < b} U {p l q h : < i < a}. Therefore, the 
maximal weight is reached when taking the whole set, and only in this case. □ 

As a consequence, a (p, q)-SCP of weight G(m) and whose parts do not exceed m 
is characterized by its greatest part only. Moreover, denoting by p a q b this greatest 
part, we have G(m) = h(a,b) : where h is the mapping defined on N 2 by 

n a+l 



h(a, b) 



q b - 1 



p" 



1 



(3) 



(4) 



q-1 p-1 

Accordingly, the definition of G may be rewritten as 

G(m) = maxfe, where P m = {(a, b) e N 2 : p a q b < m}. 

Finally observe that the greatest part of a (p, g)-SCP of weight G(m) and whose 
parts do not exceed m must be a maximal element of E n [0, m] for the divisibil- 
ity order. Otherwise, the partition could be augmented by a part, resulting in a 
partition of larger weight. The next Section is devoted to the set of these maximal 
elements. 

3. On the set Z m 

Let us denote by Z m the set of all maximal elements in E n [0, to] for the divisi- 
bility order. For convenience, we further denote by p the logarithmic ratio of q and 
p, i.e. 

logq 



P : 



\ogp 



> 1. 



Since p and q are multiplicatively independent, p is irrational. Notice that the 
elements of E may also be written as E = {p a+b P : (a, b) e N 2 }. There are exactly 
[log^TOj + 1 elements in Z m , described in the Lemma below. 

Lemma 2. Let m be a positive integer. The following characterization holds: 
p a q b £ Z rn ^==^ < b < [\og q toJ and a = [\og p m-bp\. 



4 



FILIPPO DISANTO, LAURENT IMBERT, AND FABRICE PHILIPPE 



Proof. An element p a q b of E is in Z m if, and only if, a and b are non-negative, 
p a q b < m < p a+1 q b , and p a q b < m < p a q b+1 . Since p < q, the latter condition is 
superfluous. It is easy to check that the remaining conditions are equivalent to the 
Lemma's claim. □ 



As an immediate consequence, let us note for further use that 



qZ m U {p 



h 



Z n 



pZ m if [1/p + log, mj = [log^ m\ , 

pZ m U {gL log « mJ+1 } otherwise. 



(5) 



(6) 



The elements of Z m correspond exactly to the maximal integer points below or 
on the line of equation a log p + b log q — log m = 0. An example is given in Figured] 
The corresponding values p a q b and h(a,b) are reported in Table [TJ 




Figure 2. The set Z^ for (p, q) = (2, 3), represented as all max- 
imal integer points below the line of equation x log 2 + y log 3 — 
log 750 = 0. The points along the dashed line correspond to the 
first values of the sequence £ defined in Theorem [TJ 



(a,b) 



(0, 6) (1, 5) (3, 4) (4, 3) (6, 2) (7, 1) (9, 0) 



p u, q u 

h(a, b) 



729 
1093 



486 
850 



648 
1255 



432 
850 



576 
1147 



384 
766 



512 
1023 



Table 1. The elements of Z 750 for (p,q) — (2,3), together with 
the corresponding values p a q b and h(a,b). Note that G(750) = 
/i(3,4) = 1255. 



Further define z m as the greatest integer of the form p a q b less than or equal to 
m, that is, 

z m = max Z m . 

Since qL lo s q m J g z m , we have z m — » oo when m — > oo. The next Proposition goes 
one step further. 

Proposition 2. We have the following: z m ~ m w/ien m — > oo. 
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Proof. Let z m be the smallest integer of the form p a q b greater than or equal to to. 
Thus we have z m < m < z m . By a theorem of Tijdeman [6] we know that, for to 
large enough, there exists a constant C > such that 



< 



(log z m ) c 

so that < m — z m < z m j (log z m ) c . □ 

From a greedy point of view, one might think that taking z m as the largest 
part of the (p, g)-SCP formed as in Lemma [1] yields a (p, qr)-SCP of weight G(to) ; in 
which case our asymptotics problem would be solved using the above Proposition. 
Unfortunately, this is not true. In Table [TJ we see for instance that Z750 = 729, 
obtained for (a, b) = (0,6), does not give a (p,q)-SCP of maximal weight. Instead, 
the maximal weight G(750) is /i(3, 4) = 1255. Hence, even if to is of the form 
p a q b , the first part of a (p, q)-SCP of maximal weight may be different from to. For 
instance G(729) — 1255 obtains from a unique (p, g)-SCP whose first part is 648. In 
the next Section, we study the subset Y m of Z m yielding the heaviest (p, g)-SCPs, 
i.e. those of maximal weight G(m). 



4. On the set Y m 

According to G(to) is equal to h(a,b) for certain values a, b. First notice 
that these values are not necessarily unique with respect to this property, because 
h is not necessarily one to one. As an example, in Table [TJ we observe that with 
(p,q) = (2,3) we have ft(l,5) = /i(4,3) = 850. Similarly, with (p,q) = (2,5) we 
have h(0,2) = /i(4,0) = 31. 

Let Y m be the set of all elements p a q b in EC] [0, to] such that h(a, b) = G(m). As 
already noticed, Y m is a subset of Z m . Next recall that z m = maxZ m needs not be 
in Y m . A particular link between Y m and z m however exists, given in the following 
Proposition. 

Proposition 3. For to G W , let z m = p a q b . Then 

Y m c {pV G Z m : j < b}. (7) 

Proof. Using ([3]), we have 

P -^h { a,b) = P — q ^r + U a+l -l) 

p p q - 1 p 

= p«q b -q b 



pq-p p(q-l) 



so that 



h(a, b) = — ^— (p a q b - rq b ) — , where r = E (0, lip). (8) 

p — 1 v ' q — l pq — p 

Note that < r < l/p, because pq — p > p(q —p). As a consequence, we get 

h(a,b) > h{a',b') ^ p a q b ~p a 'q b ' > r(q b - /). (9) 

If p a 'q b ' £ Z m we have z rn = p a q b > p a ' q b ' . Hence b' > b implies h(a, b) > h(a' , b'), 
which concludes the proof. □ 
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Geometrically, the above Proposition tells us that the points (a, 6) G N 2 such 
that h(a,b) = G(m) cannot be located "above" or equivalently "left" of z m . In 
particular, when z m = p a , we have G(m) = h(a, 0) = (p a+1 — l)/(f> — 1). 

In Proposition 21 we will see that the set Y m has at most two elements. For now, 
let us first focus on those elements of E that provide the heaviest (p, g)-SCP in a 
unique way, i.e. those for which Y p a q b — {p a q b }- The following Theorem shows that 
the corresponding points in N 2 form an infinite area whose boundary is a particular 
sequence, as illustrated in Figured 

Theorem 1. There exists a sequence £ — (4)beN in N such that 

Y paqb = {p a q»} <=> a > 4- (10) 
Moreover, the sequence £ is non- decreasing, unbounded, and satisfies 4 = 0. 

Proof. Let us first establish the following statements: 

(i) For all b > 0, there exists a > such that p a q b G Y paq b. 

(ii) If p a q b G Y p a qb then, for all k > 1, Y pa+kqb = {p a+k q b }. 

Let b G N. As already seen, the mapping (i,j) >-> p k qi is one-to-one. Therefore, 
we have q b — p % q l > 1 for all p l qi G Z q t \ {q b }. Choose a such that 

P a > r{q b - 1), where r = — — — as in ©. 

VI -v 

Then, for all p % q l G Z q b with j < 6, we have 

p a q b - p a+ Y > p a > r(q b - 1) > r(q b - q j ). (11) 

Using ([5]), it follows that h(a,b) > h(a + i,j); in other words p a q b G p a Zqb . -A.C - 
cording to Prop. O we also have Y paq b C {p 1 ^ G Z paq b : j < b}. Using Lemma [21 
it is immediate to check that the latter set is identical to p a Z qb , thus p a q b G Y paq b 
and (i) is proved. Now, if p a q b G Y paq b then replacing a by a + k for any k > 1 
turns (fTT|) into a strict inequality. Therefore, Y pa +k q b = {p a+k q b }, and (ii) is proved 
too. Accordingly, letting 

4 = min {a G N : Y p * qi = {p a q b }} (12) 

provides the claimed sequence I. Since Yy = {1} and 1 = p°q°, we get £o = 0. 

Let us now prove that £ is non-decreasing. Given b G N, either 4 = thus 
4+i > 4, or 4 > 1- In the latter, there exists p % q^ G Z p e b -i qb such that j ^ b and 
h{i,j) > h(£b — l,b). From ([3]), it is not difficult to see that h(i,j + 1) = q/i(i, j) + 1, 
and thus /i(i,j+l)-/i(4-l, 6+1) = q(h(i,j) - /i(4 - 1,6)). Therefore 6, (i,j+l) > 
6,(4 -1,6+ 1), so that 4 + i > 4 - 1. 

Finally suppose that £ is bounded. This would imply that there exists an integer 
a such that, for all 6 G N, Y paq b = {p a q b }. The following statement shows that this 
is impossible. 

(iii) For all a G N there exists 6 G N such that h(a + [bp\ , 0) > h(a, 6). 

Indeed, by Lemma El we know that p a +\- b P\ e Z paqb . Fix a G N, choose 6 G N*, 
and set a 1 = a + L6pJ, where p = logq/logp. According to (|9|), and denoting by 
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{bp} — bp — \bp\ the fractional part of bp, we have 

h(a, b) < h(a', 0) e> p a q b - p a ' < r{q b - 1), 

«*■ p LbpJ > g h (1 - r/p Q (1 - l/g 6 )) , 
«• {&/>} < - log p (1 - r/p Q (1 - 1/g 6 )) . (13) 
Now observe that sequence (<fi a ,b)b&t defined by 

^ = -^(l-^(l-^)) (14) 

is increasing, with ^> aj o = 0. Since p is irrational, we know that (bp)b&i is equidis- 
tributed modulo 1. Thus, there exists b > such that {bp} < 4> a .i, hence {bp} < 
tfia.b, which using fT3")) concludes the proof. □ 

Let us anticipate a result of the next section, implying that the sequence I is 
completely known as soon as we can compute its jump indices, that is, the values 
b > for which lb > lb-\- Indeed, we shall establish with statement (|25|) that if b 
is a jump index of I then lb = L a (^)J 1 where 

^ lQ g P \ z p g { bp} + lQ gp f ■ ( 15 ) 

Computing the jump indices of I may be done recursively as shown in the next 
Corollary. Referring to the sequence cf> defined in let the mapping j3 be defined 
on N by 

(3(a) = min {j S N : j > 0, {jp} < OJ -}. (16) 

Corollary 1. T/ie increasing sequence (bk)keN of the jump indices of £ satisfies 

6 = /3(0), 6 fc+ i=/3(4J. 

Proof. Given any a S N, we know that Y p a = {p a } by Prop. [3J Moreover, letting 
b be incremented by 1 from iteratively, it follows from ([5]) and © that Y paq b — 
{p a q b } as long as h(a,b) > h(a + , 0). As it is shown in part (hi) of the proof 
of Theorem [TJ the latter inequality is equivalent to {bp} < <f> a .b- Accordingly, 
Y paq b = {p a q b } if b < /3(a). Hence £p/ a )-i < a < £pt a \, so that /3(a) is a jump index 
for £. The result follows immediately since I is non-decreasing. □ 

As claimed before, we next show that Y m has at most 2 elements. This might 
be established by directly using ^ and studying the diophantine equation 

p a q b - p c = r(q b - 1), where r = -iZil. 

pq-p 

Unfortunately, the latter is seemingly not easy to cope with, whereas Theorem Q] 
proves handy. 

Proposition 4. For all m S N, the set Y m has either one or two elements. 

Proof. Assume #Y m > 2 and denote by p a q b its greatest element. Since Y m = Y paq b, 
we have a < lb by definition of sequence I in Theorem [TJ To be more precise, 
statement (ii) in the proof of this theorem even tells us that a — lb — 1 . Now let 
p°q d be the second greatest element in Y m . According to Prop. Owe must have 
d < b, thus c > a by Lemma [21 Since £ is non-decreasing, it follows that Id < 4- 
Since a = 4 - 1 we get c > £<j, which means that Y pCq d = {p c q d } from Theorem [TJ 
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Therefore, there cannot be a third element in Y m as it would also be an element of 

Ypc q d. □ 

5. Asymptotic behavior of G 

Our goal in this Section is to prove that G(m) is equivalent to mpj (p — 1) as to 
tends to infinity, independently of q. As a simple first step, let us exhibit a sharp 
upper bound for G. 

Lemma 3. For all to £ N and all n £ Y m , we have G(m) < npj(j> — 1). In 
particular, 

G(m) p 
hm sup = . 

m^oo m p - 1 

Proof. Let n — p a q b £ Y m . According to ([5]) we have 

h(a, b) ~ =-( + < 0, where r £ (0, lip). (17) 

P-1 \p-l q-lj 

Hence G{m) = h(a.b) < np/(p — 1). Since n < to, it follows that G(m)/m < 
pj (p— 1). To conclude, observe that for to = p a we have G(p a ) = (p a+1 — l)/(p— 1) 
from PropO Therefore, \im a ^ OQ G(p a )/p a =p/{p— 1)- □ 

Let us now define a mapping y as follows: For all to, let y m denote the smallest 
integer of the form p a q b such that G(m) = h(a, b), that is 

y m = m\nY m . (18) 

According to Proposition [31 y m is also the element of Y m with the smallest 
exponent in q. We shall next give a characterization of y m using the sequence 
I defined in Theorem [1] Recall that this sequence is defined by lb — min{a e 
N : Y p a qb = {p a q b }} and satisfies (fTTJl) . Since I is non-decreasing, the sequence 
(p eb q b )befii is increasing. We may thus define, for all to € N, 

m e = max{6 £ N : p lh q b < to}. (19) 

Theorem 2. For all m £ N, we ftawe 

y m = max{p a 9 h £ Z ro : b < m e }. (20) 

Moreover, let a — [log p to — TO^pJ and to = \ m/p a \ . Then y m — p a Zm- 

Proof. Let y m — p % q>. Since y m = minY^, we have Y Vm = {y m } thus i > £j 
using (fT0|) . Suppose j > mi, then > p ij q^ , and thus > to from (fT9|) . 
which contradicts the fact that y m < to. Therefore, j < m£. 

Now consider any p a q b £ Z m such that 6 < to^. Since £ is non-decreasing we 
have i mt > ^6. By Lemma HJ there exists A: £ N such that p k q me e Z m , and we 
have fc > i mi because p lm eq m t < m . Since p a g b £ Z m , condition b < me implies 
that a > k by Lemma [21 again, so that a > £ me > h- Therefore, Y paqb = {p a q b }. 
Supposing p a q b > p i q l would then imply that h(i,j) < h(a,b), contradicting the 
definition of y m . Thus p a q b < p l q 3 , and ([2TJ]) is established. 

Accordingly, j = [log p TO — ip\ < a, so that y m /p a is an element of E. Therefore, 
y m /p a 5: m/p a implies y m /p a < \_ m /p a \ = m i which in turn implies \m/p a \ < z r7l . 
We thus get 

Vm < p a Zm < p a m < TO. 
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Let Zf- n — p a q b . To conclude the proof by using (j2"0|) again, it suffices to show that 
b < me. Since p d q me G Z m , we have p"q me <m< p "q™t+ 1 , so that q mi < m < 
q m i+ l , Thus \}og q m\ = mi, whence b < mi by Lemma [5] □ 

Comparing with Proposition [3J characterization (|2U|) of y m no more depends on 
z m . Moreover, it provides a first improvement of Lemma [JJ 



Corollary 2. For all m G N, we Ziawe 

G(m) p 



as m — > oo. (21) 



Z/m P - 1 

Proof. For m G N, let y m = p a,n q bm . According to (|17p we have 

G(m) p k m q-p 

— where k„, = 



Vm P-I P a ™' (p-l)(?-l) q b ™(q-l) 



Observe that k rn G y (p-i)(g-i) ' jAl * s UIU f° rrru y bounded. To conclude the proof, 
we need to prove that a m -4 oo as m -) oo. According to Th.[21 we have b m < mi, 
thus a m > £ me . Since mi goes to infinity with m, and since I is non-decreasing and 
unbounded by Th. [1] a m goes to infinity with m too. Hence the claim. □ 

According to the latter result and Proposition^ the final task consists in showing 
that y m ~ z m . This is done next, so that our main claim is established. 

Theorem 3. For all m G N, we have 

G(m) p 



to p — 1 

Proof. Assume y m ^ z m . According to Th. [21 the elements in Z m that exceed y m 
are of the form p l q 3 with mi < j < Ll°g g m J ■ Let us sort these elements together 
with y m in an increasing sequence (n,i) i£ [ ^, so that no = Um and 77 at = z m . 
Notice that the elements of this sequence are consecutive elements of E for the 
usual order. As soon as m is large enough, we know by Tijdeman's result already 
mentioned [5] that 77, + i — rii < n,i/(\ogni) for an explicitely computable constant 
C > 0. Therefore, 

N-l JV-1 

(log 77^ ^ 

Accordingly, for all m G N, 



< Z m - y m < ([l0g ? 77T.J - 777^) ^mjC ' ( 22 ) 

At this point, what remains to be proved is that \\og (j m\ — mi grows asymp- 
totically slower that (logy) so that z m ~ y m , and to conclude using (|21|) and 

Lemma [5] Recall that 777^ is defined as the largest value b such that p lh q h < m. 
Therefore, we have 

or equivalently, using p = log 5/ logp, 

— - + 777^ < log 777 < - \-mi + l, (23) 

P P 
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SO that 

Llog„mJ-m £ <^±i + l. (24) 

P 

It thus remains to evaluate the terms in sequence £. For that purpose, we first 
give an explicit formula for £, valid at the jumps of I. We claim that for all 6 6 N*, 

4>4_! ^=K j!4 (& -p^) J- (25) 

Indeed, assume that 4 > 4-1- Then, for all p l q> € Z p e b _ lqb _ 1 , we know from @ 
that 

p^q b - 1 -p l q> > r{q b ~ l -q>). 
Multiplying both sides of the above inequality by q yields, for all p l q J £ qZ p i b _ lqb _ 1 , 

V tb -^q h -pV > r{q b - qi). 

Now, 4 > 4—1 implies that there exists an element p l q J G Z p e b _ lqb for which the 
latter inequation does not hold. By Lemma [Hand the definition of Z qm in this 
element must be p e >>~i + [ b Pl ; so that 

p^-i^-pLM) <r(q b -l). 

In fact, note that this inequality does not only hold for p lb - x ; by definition of £, it 
remains valid for jr 1+1 , . . . ,p ff > — 1. Accordingly, we get 

p^-^qb _ p L6pJ ) < r ( g b _ i) < p h ( q b _ p [b P \ ^ 
which proves claim (|25[) . 

It follows from ([23]) that, for any 6 such that 4 > 4-ij 

Using another result of Tijdeman (see Th. 1 in j7]), we know that, as soon as 
pY b p\ > 3 ; there exists another explicit constant C > 1 such that 

g 6_pLhpJ> 1 (27) 

" (log P L6pJ) c 

Therefore, since g b = p bp , and (g7J) imply that 

„bn n „ \bp\\C' fii 

4 < log p q 1 g j ? > = {>} + logpLM < l + -i— log log <z b . (28) 

^ pL 6 PJ logp logp 

Putting all this together, we get the claimed result. Indeed, let b be the smallest 
index such that £ me +i = 4- Since 4 > ib-i, we have 

c c c 

im e +l = 4 < 1 + 1 log log q b < 1 + log log q me+1 < 1 + log log qm. 

logp logp logp 

(29) 

Using ([M]) we get 

1 C i t\ 

I log ml - m,£ < 1 H h log log = o ( (log m/p) c ) , (30) 

P logg V J 

which implies, using ([22]) . that y m ~ z m and concludes the proof. □ 
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Note that the above proof mainly relies on the fact that the sequence i is non- 
decreasing and that, due to the lower bound in (|27j) essentially, it grows very slowly. 
The theorem of Tijdeman that provides this lower bound hinges on a result of 
Fel'dman about linear forms in logarithms. More recent results of Laurent et alii [5] 
about such forms in two logarithms allow one to make precise the value of the 
effective constant C in (f^Tj) . Nevertheless, the algorithm presented in the next 
Section does not make use of ([3"0"|) in order to compute mi in logarithmic time. 

6. Computing y m and G(m) in log log m time 

Making use of the mapping h, computing G(m) is straightforward as soon as 
an element of Y m is known, in particular y m . According to Theorem [21 y m /p a is 
the greatest element of Z fh , where a and fh depend on wig. Once the latter value 
is known, computing the greatest element in Z„ can be done in logarithmic time 
with an algorithm explained in [9 . 

According to Corollary [TJ mi may be effectively computed if an efficient way of 
computing the mapping (3 is found. We shall see that the (3(a) 's are denominators 
of convergents of p, which already improves the computation time. But we shall 
also see that the latter property implies that the relation lb = [ct(b)\ , see ([15]), also 
holds for all denominators of both even convergents of p and their mediants. This 
allows computing me even more quickly by only using a simplified version of the 
mapping a and binary search. 

Note that the convergents of p are also of use for the algorithm in [5]. Let us 
first recall some well known facts about them (see, e.g., [10] or [5]). Let [ao,ai, ...] 
be the regular continued fraction converging to p, and denote by hi/ki the i th 
principal convergent of p (i > 0). The sequences (/i2i/fc2i) and (/i2i+i / 'k2i+i) are 
adjacent, converge to p, and {ki + ki+i)^ 1 < \kip — hi\ < h l+ \ . Given i > 0, 
the intermediate convergents of hi/ki, sometimes referred to as mediants, are the 
rational numbers hij/kij, for < j < a,i + 2, given by 

h,j = K + jh l+1 , k hj = k t + jk l+1 . (31) 

Let us denote by (H n / K n ) n( z^ the increasing sequence of all convergents of p of 
even index together with their intermediate convergents. It is known ([11], Th. 2) 
that (H n /K n ) is the best lower approximating sequence of p, that is, its terms are 
characterized by the following property. For each n6N and integers h, fc, 

H n h 

-^<-<p k>K n . (32) 

Two immediate consequences are needed here. The first one is that, while the 
sequence (K n ) increases to infinity, the sequence {{K n p}) decreases to 0. Indeed, 
property ([3"2"| implies that H n = \K n p\, so that (|3"T]) implies, for < j < a.i + 2, 

fej+ip} - {hi^p} = k 2 i+iP - h 2l+ i G {-l/k 2i+2 , 0). (33) 

The second one rephrases the sufficient condition in ([32]) : If {bp} < {jp} holds for 
all integers j such that < j < b, then 6 is a term of (K n ). Indeed, let h, k be such 
that 

[bp\ h 
— < k <P - 
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Since h < kp, the above inequalities still hold for h — [kp\ , which implies {kp} < 
■|{6p}. Supposing k < b yields {kp} < {bp}, a contradiction, so that \ bp\/b satisfies 

(E2D- 

Now our previous claims may be established. 

Theorem 4. For all a G N, (3(a) is a term of the sequence (K n ) n eN- Moreover, 
for each K G (K n ) neN , t K = \a(K)\ . 

Proof. Since the sequence (4> a ,i) increases from with i and the sequence ({Kip}) 
decreases to 0, there is a unique n such that {K n p} < <p a ,K n and {Kip} > 4> a ,Ki 
for all i < n, if any. In particular, K n > (3(a). We next establish that K n = (3(a). 
This clear if n = since Kq = 1 and /3(a) > 1, so assume n > 1 in the sequel. 

Let b = (3(a) for short, and suppose b < K n . Let K = K n -\ for short again, 
and let c = min{j > : {jp} < {Kp}}. For alH < c we have {ip} > {Kp} > {cp}, 
so that c S (Ki), which implies c = K n since (Ki) decreases. Therefore, b < K n 
forces {bp} > {Kp}, so that 

K.b > {bp} > {Kp} > c(> atK . (34) 

Since (f) at i increases with i, we get K < b < K n . Property ([32]) implies that 
[bp\/b< \K n p\/K n . Since we cannot have [Kp\/K < [bpj/b < [K n p\/K n f[TTj. 
(ii) of Lem. 1), it follows that \bp\/b < [Kp\/K, that is, {bp}/b > {Kp}/K. Thus 

fh\ 1 . b - K <r, i ^ Kk v a(l-q- K ) 

{bp} {Kp} > —{Kp} > — > KlQgp ■ 

Nevertheless, 



^a,h - ^a,K < <^a,oo - <f>a,K = log„ I 1 + " Q K ] < , , 

s 1 — " ' (1 — oe)q K iogp 



According to (|34|) we should thus have 

(i-q- K ) , 

K (l-a)q K ' 
which would imply 

q K - 1 1 p(q-l) 

K 1 — a q(p — 1) 

Therefore, K n — (3(a) as claimed, which proves the first assertion of the Theorem. 

Next turn to the second one, and let K € (K n ). One the one hand, Ik < [o;(-K')J 
holds. Indeed, there is a unique jump index K* of I such that Ik = £k* , and K* < 
K because I is non-decreasing. According to the first assertion of the Theorem, 
K* G (K n ). Since (K n ) increases and ({K n p}) decreases, (a(K n )) is increasing, 
thus (\_a(K n )\) is non-decreasing. Therefore, Ik — = [a(K*)\ < [a(K)\. On 
the other hand, we also have [a(-?0J < £k- Indeed, letting a = [a(K)\ for short, 
we have a < a(K), that is, 

na< ( q ~ P )(q K -l) 
P ~ (q-l)(q K -pL-M)' 

Recalling ((§]), this also reads 

(p a - 1 q K -p a - 1+lKpi ) < r(q K - 1). 
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Since £k > 0, we may assume a > 1. Letting a' = a + \Kp\, the above inequality 
means that 

h(a-l,K) < h{a! - 1,0). 

Finally notice that p a ' _1 < p a ~ 1 q K . Therefore, Y p a-i q K {p a ~ 1 q K }, so that 
a - 1 < £ K by ([TP]). Whence [a (if) J = a < £ K as claimed. □ 

We may thus use the following algorithm in order to compute me given m, p, 
and q. 

1 . Compute the denominators k t of the primary convergents of p until finding 
the largest one with even index, say k* = fcaij that satisfies p k * q k < to. 

2. Compute the largest j such that K* = kiij satisfies p ll< * q K < to. 

3. Return mi = \\og q m — Ik* / p\- 

Let us finally evaluate the number of operations required for these computations. 
Since fej > 2 2 , Step 1 requires computing O(loglogm) values. Each of them ob- 
tains in O(l) time using any standard formula for generating the dj's and the k^s. 
Checking the stopping conditions in Steps 1 and 2 requires computing = a(fc), 
but this may be done in 0(1) time without computing q k . Indeed, let 

a(b) = \og p iZ£ _ lo gp (l -p-i^) = a(b) - log p (l - q- b ). 

Since a(b) — a(b) lies in the interval (^-, - q irzi), ^ or eacn K ^ (^n) we either have 
Ik = a{K) or i 1 ^ = &{K) — 1. Therefore, any if satisfying p a+< - K ^q K > m also 
satisfies p lK q K > to. So a is only needed to compute the true value of I two times, 
in order to check and Ik*- This is done in O(loglogTO). Finally, there might 
be a large number of j-values to be coped with in Step 2, but this can also be done 
in O(loglogTO) time using binary search. 
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